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Abstract 

In this paper, we propose a reformulation of the bootstrap Multi-Scale Analysis (BMSA), 
developed earlier by Germinet and Klein, as a single scaling algorithm, to make explicit the 
fact that BMSA technique implies an asymptotically exponential decay of eigenfunctions 
(EFs) and of EF correlators (EFCs), in the lattice Anderson models with diagonal disorder, 
viz. with an IID random potential. We also show that the exponential scaling limit of the 
EFs and EFCs holds true for a class of marginal distributions of the random potential with 
regularity lower than Holder continuity of any positive order. 

The version info: compared to v.l of 9 March 2015, in the v.2 (21.07.2015) the scale-free 
threshold has been slightly improved: 841 d = 29 2d , as in [28], is replaced by 23 2d , owing to an 
optimization of the geometrical argument in the proof of Lemma 1; the present text contains ad¬ 
ditional figures; some intermediate technical parameters have been eliminated; all the remaining 
parameters and the relations between them are listed in the tables (4.16)-(4.18) [(4.13)-(4.14) 
in v.3]; the proof of the key Lemma 6 is rendered more detailed. 

v.2 —> v.3: Fig.6 has been added, and some numerical parameters has been rectified in Fig.5. 


1 Introduction and motivation 

We consider Anderson models with diagonal disorder in an integer lattice Z d , d > 1. Such mod¬ 
els have been extensively studied over the last thirty years; the two principal tools of the modern 
rigorous Anderson localization theory are the Multi-Scale Analysis (MSA) and the Fractional 
Moment Method (FMM). In the framework of lattice systems (and more generally, systems on 
graphs with sub-exponential growth of balls) the MSA proved to be more flexible; in particular, 
it is less exigent to the regularity properties of the probability distribution generating the local 
disorder - in the simplest case, the single-site marginal distribution of the IID (independent and 
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identically distributed) values of the external random potential. On the other hand, a consider¬ 
able advantage of the FMM in the same class of models is to provide exponential decay bounds 
for the (averaged) eigenfunction correlators (EFCs), under the condition of Holder continuity of 
the single-site marginal distribution. By comparison, the original MSA scheme by Frohlich et 
al. [23], reformulated by von Dreifus and Klein [19], proved only a power-law decay of the key 
probabilistic estimates in finite volumes. When the MSA was adapted to the proofs of strong 
dynamical localization (cf. [25, 20, 26]), this resulted in power-law decay of EFCs. 

Germinet and Klein [26] significantly narrowed the gap between the EFC decay bounds 
provided by the MSA and FMM. Specifically, using the bootstrap MSA, involving several in¬ 
terconnected scaling analyses, they proved sub-exponential decay bounds with rate L H > e _L<s 
for any 5 G (0,1). Their result is a finite-volume criterion of exponential spectral and sub¬ 
exponential dynamical localization: if the latter occur for a random Hamiltonian of one of the 
forms considered in [26], then their scale-free condition is fulfilled for L 0 large enough. 

Recently Klein and Nguyen [28, 29] have adapted the BMSA to the multi-particle Anderson 
Hamiltonians. 

In theoretical physics, the celebrated scaling theory, put forward by the ”Gang of Four" 
(Abrahams, Anderson, Licciardello and Ramakrishnan, [2]) and further developing the Ander¬ 
son localization theory [1], predicted - under certain assumptions including also those sufficient 
for the MSA or FMM to apply - that the functionals F L related to the quantum transport, first 
of all the conductance, for systems of large size L, should admit a limiting behaviour in the 
double-logarithmic coordinate system. While the existence of a.c. spectrum for systems on a 
periodic lattice or in a Euclidean space remains an intriguing challenge for the mathematicians, 
we show that in the parameter zone(s) where various forms of localization can be established 
with the help of existing techniques, the rate of decay F(L) of eigenfunction correlators (EFCs) 
at large distances L admits the limit 
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Below we will call such a behavior exponential scaling limit (ESL). Formally speaking, we 
obtain, as usual, only upper bounds, but the example of one-dimensional systems shows that 
decay faster than exponential should not be expected. 

The main goal of the present paper is a transformation of the Germinet-Klein multi-stage 
bootstrap MSA procedure from Ref. [26] into a single scaling algorithm, replacing several inter¬ 
connected scaling analyses in the bootstrap method and establishing the ESL in the traditional 
Anderson model. 

The motivation for the present work came from an observation, made in Refs. [11] (cf. [11, 
Theorem 6]), [13] (see Theorem 8 in [13] and discussion after its proof), and some earlier works, 
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that already in the von Dreifus-Klein method from Ref. [19] there were some unexploited 
resources, giving rise to “self-improving” estimates in the course of the induction on the length 
scales L k , k > 0, following the recursion L k = [L k _ L \ ~ Lq , a > 1. Specifically, it was 
observed that the k -th induction step actually produces more decay of the GFs than required for 
merely reproducing the desired decay rate at the step k + 1, and that this excess can be put in 
a feedback loop, improving the master parameters of the scaling scheme. The net result is the 
decay of the GFs (and ultimately, EFCs) faster than any power law 1 , viz. L (->■ e ~ ahil+c l, with 
a, c > 0. 

The benefits of such a feedback-based self-enhancement of the master scaling parameters 
become much greater when the scales grow multiplicatively, as in the first stage of the BMSA: 
L k = YL k _i = Y k L 0 , with Y > 2. A fairly simple calculation shows that essentially the same 

k 

feedback loop as the one used in [11, 16, 13] for the scales L k ~ Lq , k > 0, gives rise in this 
case to a fractional-exponential decay L (->■ e~ LS , with some <5 > 0. 

Acting in the spirit of the bootstrap MSA, we implement a technically more complex scaling 
procedure than the above mentioned “simple feedback scaling”, aiming to render more explicit 
and constructive the statement of the BMSA (cf. [26]) that any (viz. arbitrarily close to 1) value 
of the exponent 5 in the above formulae can be achieved for L large enough. To this end, we 
replace the first two stages of the BMSA (with fixed parameters) by an adaptive feedback scaling 
algorithm. The latter makes the multiplicative growth factor Y, figuring in the scaling relation 
L k = YL k _i, scale-dependent: Y k = y(k. L k ). In fact, the BMSA scheme includes another 
important geometrical parameter - an integer S k G [1, Y k ); see Section 3. 

However, the ’’simple feedback scaling” - with Y k and S k fixed - may still be required during 
an initial “boost” stage, where the effects of localization are almost imperceptible, particularly 
in the probabilistic estimates. Since the scales grow with k (viz. L k = Y k L 0 ), writing formally 
Y k = L T k results in a finite, initial sub-sequence {ti, ..., r^_i}, with some Si depending upon the 
model parameters, which is actually decreasing. (As such, the values {ti, ..., t$_i} are simply 
unused.) It is only later, for k > that we fix r k = r > 0, thus effectively switching to the 

Cl_i -j-\k 

super-exponential growth L k ~ CL^ .Of course, depending on the reader’s personal point 
of view, the presence of this switching point may be considered as a form of the Germinet-Klein 
multi-stage technique. 

Taking into account the abundance of various technical parameters in our scheme, we keep 
r k fixed for the rest of the scaling procedure. However, the algorithm’s efficiency can be further 
improved by making r k also ^-dependent (and growing). This may prove useful in a numerical 
implementation of the adaptive scaling algorithm, as well as in specific models (including the 

'This result holds true under a very weak regularity of the random potential, just barely stronger than the 
conventional log-Holder continuity of the marginal distribution. See Assumption (W3) (Eqn. (2.15)) in [13]. 
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multi-particle models with slowly decaying interaction). We show that the ’’gap” between the 
genuine exponential decay (i.e., the value 5 = 1) and the exponent 5/, achieved at the A'-th 
step, decays at least exponentially fast in k. In a way, this provides a rigorous complement 
to the predictions of the physical scaling theory on the convergence to the ESL, at least in the 
parameter zone(s) where localization can be proved with the existing scaling methods. 

The core of the renormalization group-style analysis, the proof of the key Lemma 6, is not 
very long, but tedious, and the number of intermediate scaling exponents and other parameters 
is quite disturbing, so the reader may naturally wonder if the entire inductive algorithm could 
be made significantly simpler and more transparent. The answer is affirmative, but this would 
require a certain concession: as show our preliminary calculations, it suffices to start with the 
initial probability p 0 (cf. (4.11)) which is not scale-free but polynomial in L 0 , i.e., po < 
with a value of fro > 0 which is, actually, not excessively large. This would eliminate the need 
for the additional ’’boost” phase in the scaling algorithm, where the decay of the probabilities 
of unwanted events is barely perceptible (and virtually useless for applications to physically 
realistic 2 * * models), and also render the whole procedure more straightforward. In the author’s 
opinion, this would be indeed only a small concession, for two reasons: 

• the main role of the scale-free Germinet-Klein’s criterion is to provide a mechanism of 
transforming weak probabilistic bounds into much stronger ones, and this task is entrusted 
to the adaptive feedback that we propose; 

• more importantly, the actual verification of the initial conditions for the onset of An¬ 
derson localization is usually based on one of the two methods: using the Lifshitz tails 
phenomenon (providing de facto very strong, fractional-exponential bounds), or deriving 
the initial bounds from the strong disorder assumption, which also provides strong scale- 
dependent probabilistic bounds. 

Still, we keep in the present paper the scale-free localization criterion to show that the 
Germinet-Klein BMSA extends to the result on exponential scaling limit without any additional 
hypotheses. 

Speaking of the consecutive phases (analyses) of the BMSA, it is to be pointed out that 
we do not perform the last stage where a genuine exponential decay of the Green functions is 
established in cubes of size Lp- with probability ~ 1 — e~ k , where 5k = 5 is made arbitrarily 
close to 1 by the results of [26]; one would expect 1 in the framework of the present paper. 
We do not analyze the behaviour of such probabilities related to the exponential decay of the 

2 To have an idea of typical situations where Anderson localization theory proves valuable in modern technology, 

recall that a typical CPU is today a square film with about 10 7 :- 10 s atoms alongside (~ 1cm) and just a few atoms 

(~ 20 nm) across it. 


4 



GFs in finite volume. As was already said, this paper focuses on the exponential scaling limit 
- for the Green functions, eigenfunctions and eigenfunction correlators. The actual road map 
is as follows: GFs -w EFCs EFs, so the decay rate of the EFs is shaped by that of the 
EFCs. Naturally, one can switch at any moment from the analysis of the ’’almost exponential” 
decay to that of the exponential one, by simply following the Germinet-Klein approach, but our 
main goal is the construction of a single algorithm which takes care of all exponents 5 close to 
1. Undoubtedly, there can be various further developments of the BMSA technology from [26]. 

Finally, we show that the proposed adaptive feedback scaling technique allows for a lower 
regularity of the marginal distribution of the IID random potential than Holder continuity of any 
positive order. In the realm of the FMM proofs of localization, it is well-known that absolute 
continuity of the marginal distribution can be safely and easily relaxed to Holder continuity 
of any positive order f3 (cf. [3]) Yet, the MSA in general is renowned for its higher tolerance 
to a lower regularity of the probability distribution of the disorder. So, while the question on 
the lowest regularity compatible with the FMM approach to the exponential strong dynamical 
localization remains open, our results evidence that Holder continuity is not required for the 
exponential scaling limit of the EF correlators. 

As was said, strong dynamical localization at some fractional-exponential rate 5 6 (0,1) 
actually follows from the initial, weak hypotheses through a simpler scaling procedure, under 
the assumption of Holder continuity of the marginal PDF of the random potential. 

We would like to make one last comment, hopefully providing an answer to the readers, 
familiar with the FMM approach, who may wonder: does one really need an alternative method, 
such as (a version of) the MSA, to prove ’’almost” exponential localization in the models with 
sufficiently regular probability distribution of the random potential, where the FMM proves a 
genuine exponential decay of EFs and of EFCs? Well, there are several elements that have to be 
taken into account. 

• Firstly, the MSA or its variants have been successfully applied to the models with de¬ 
terministic disorder where the dependence between the values of the potential at distinct 
points is so strong that no decoupling inequalities or similar techniques, which have been 
the cornerstone of the Aizenman-Molchanov method since the pioneering work [3], seem 
to apply. 

• Secondly, the Fractional Moments Method ceases to be a ’’mono-scale” technology in the 
area of multi-particle Anderson localization. Indeed, both the paper by Aizenman and 
Warzel [7] and a recent work by Fauser and Warzel [24] are based on a scale induction for 
the fractional moments of the Green functions. This is unrelated to the regularity of the 
marginal distribution. 
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• Thirdly, in the iV-particle Anderson models with interaction decaying slower than expo¬ 
nentially, the only known proofs of genuine exponential decay of the EFs (cf. [17, 14]) 
employ an important particularity of the MSA: the analyses of the EFs and of their cor¬ 
relators can be carried out independently, so one can prove exponential decay of the EFs 
even in some situations where the EFCs decay sub-exponentially. In the framework of the 
FMM, the analysis of the EFs is subordinate to that of the EFCs. Again, this issue arises 
even in the most regular models of the disorder. 

Summarizing, mathematical methods are not always ordered on a linear scale, and new chal¬ 
lenging models may test and contest the relations between them. 

1.1 The model 

We focus on the case where the configuration space of a quantum particle in an external random 
potential is the lattice d > 1, and consider the random Hamiltonian H{u) of the form 

(H^)(x)= (V’Oc) -V’(y)) +V(x-,u)^(x), (1.1) 

\y-x\i=l 

where V : 7L d x Q —* M is an IID random field relative to some probability space ((],£, P), 
and |x|i := |xi| + • • • + \x d \ for x = (aq,..., x d ). Until Section 7, we assume that marginal 
probability distribution function (PDF) F v , of the random field V, 

F v {t) := ¥ {V(0]uj) < t} , f el, 

is Holder-continuous of some order (3 E (0,1). In Section 7 we show that the assumption of 
Holder-continuity can be slightly relaxed (cf. Eqn. (7.1)). 

The second-order lattice Faplacian in (1.1) can be easily replaced by any (self-adjoint) finite- 
difference Hamiltonian of arbitrary finite order, without any significant modification of our algo¬ 
rithm. Indeed, we replace the form of the Geometric Resolvent Inequality most often employed 
in the MSA of lattice models, with its variant traditional for the MSA in continuous systems (in 
W l ). It is based on a simple commutator relation, so that the range (order) of a finite-difference 
kinetic energy operator becomes irrelevant (and unused in the intermediate calculations), pro¬ 
vided the initial length scale L 0 is large enough. For clarity, we work only with the standard 
lattice Faplacian. 

1.2 The main assumption 

The principal assumption on the parameters of the model at hand is the probabilistic inequality 
(4.4) which we formulate in Sect. 4.1, for it requires several definitions given there. 
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1.3 Structure of the paper 

• The principal objects and notations are introduced in Section 2. 

• In Section 3, we present the main analytic tool of the scaling analysis - the Geometric 
Resolvent Inequality (GRI), and formulate the main result of the paper, Theorem 3.1, 
providing a scale-free criterion for the exponential scaling limit of the Green functions. 
The exposition is closer to the form of the GRI used in the continuous systems than to the 
one traditionally used in the lattice models, starting from the pioneering papers [22, 23, 
18, 33, 19]. This is required for the geometrical optimizations a la Germinet-Klein [26]. 

• The core of the paper is Section 4, and the staple there is Lemma 6 concluding the fixed- 
energy MSA (FEMSA). 

• The derivation of the exponential scaling limit from the results of Section 4 is given in 
Section 5. 

• Section 6 is devoted to a ’’soft” derivation of the variable-energy MSA (VEMSA) and of 
strong dynamical localization from FEMSA carried out in Section 4. 

• In Section 7, we show that the assumption of Holder-continuity of the marginal probability 
distribution of the random potential an be relaxed. To the best of the author’s knowledge, 
this result is new. 

• A number of definitions of various technical parameters and relations between them are 
listed in Sect. 4.2, in the tables (4. 13)— (4.14). The reader may find it helpful to have them 
printed on a separate page when checking the proofs. 

In theoretical physics, a sufficiently fast decay of the Green functions away from the diag¬ 
onal is usually considered as one of equivalent signatures of Anderson localization. Speaking 
mathematically, this is a higher-dimensional analog of positivity of Lyapunov exponents in one¬ 
dimensional (or quasi-one-dimensional) systems. While it is known that, in general, this analog 
does not necessarily imply spectral localization, first, it has been shown long ago by Martinelli 
and Scoppola [30] that it rules out a.c. spectrum with probability one, and secondly, it has been 
observed that the s.c. spectrum occurs in systems with some strong “degeneracies” in the prob¬ 
ability distribution of the ergodic (not necessarily IID or weakly correlated) potential. Under 
reasonable assumptions on regularity of the IID random potential, fast decay of the GFs implies 
indeed spectral and strong dynamical localization, and the role of Section 6 is to summarize 
the progress achieved in this direction, and to show in a fairly simple way that the fixed-energy 
analysis is the heart of the localization analysis of the conventional lattice Anderson model. 
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The continuous systems are not considered in the present paper, since the analysis of un¬ 
bounded (differential) random operators would require an additional technical discussion per¬ 
taining to the domains, self-adjointness, etc. But as was already said, here we focus mainly on 
the scaling algorithm that could be applied, essentially in the same way, both to the discrete and 
continuous systems. 

2 Basic geometric objects and notations 

Following essentially Ref. [26] (where the Anderson-type models in a continuous space W 1 
were considered), we work with a hierarchical collection of lattice cubes, with specific centers 
and positive integer side lengths L k . For our purposes, it is more convenient to start with the 
cardinalities of the cubes and those of their one-dimensional projections: we fix odd positive 
integers Y > 1, 4 and set 


L k = Y k • 34 = 3 ■ Y k l o =: 34- 


(At some moment in the proofs, the scaling factor Y becomes variable: Y = Y k .) Next, we 
consider the lattice cubes with coordinate projections of cardinality L k : 



Since L k = 3 Y k £ 0 is odd, the upper bound in the above definition of the cube B Lk ( x ) could have 
been replaced with (L k — l)/2, resulting in the same lattice subset. However, having in mind 


the canonical embedding 7L d W 1 , the above definition looks more natural when transformed 
as follows: with y e 7L d ^ M d , 



so that the ’’fictitious” radius of the ball is precisely L k /2. 

Sometimes it is more convenient to refer to the spherical layers and balls relative to the 
max-distance, with a clearly identified integer radius : 


C r {u ) = {i6Z' 1 : \x — u\ — r } , 

A r (u) = { x G : \x — m| < r } = B 2r +i(w) • 


( 2 . 1 ) 


Notice that one has B L k (u) = A & fc -i (u). 


2 


The cube B Lk {u) is partitioned into 3 d adjacent cubes called /.'-cells, 


Cfc(c) := B 4 (c) = Ay-i(c) 

2 


( 2 . 2 ) 


(recall: 34 = 4,) with centers c in the sub-lattice (3Z) d . 
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Figure 1: Cubes and cells. Here d = 2, Lq = 9. 3 2 — 1 cells forming the shell of a cube Bg(-) is shown in gray 
color, and the core in dark gray. 


• The central cell C k {u) of a cube B Lk ( u) will be called the core of B Lk (u); 

• the complementary annulus, formed by the remaining 3 d -1 cells of B Lk (u), will be called 
the shell of B Lk (u). 

Given any length scale L k = Y k L 0 , we shall always work with the family of L/.-cubes 
whose cells form the uniquely defined partition of Z d including the cube centered at the origin, 
B/ (: (0); these cores, as well as their centers, will be called admissible at the scale L k . The 
centers of the admissible £ fc -cores form a sub-lattice of Z d denoted by C k . Sometimes we use 
notation (c, c'), meaning that c, d G C k are two nearest neighbors (in C k ) relative to the max- 
distance: | c — d | = i k . By a slight abuse of notations, we will write, e.g., X^( C c')ec fc i nstea d of 
£ (cc , )eCfcxCfc Each point c G C k has 3 d — 1 nearest neighbors (within C k ). 

See Fig. 1 where 

• an admissible square of size L = 9 (thus with 9 vertices along each side) is shown in gray 
color; it is partitioned into 3 d = 3 2 congruent cells separated visually by thin white lines; 
the admissibility means that the periodic sublattice of the cell centers (large black dots) 
includes the origin 0 G Z d ; 

• the core, i.e., the central cell, is shown in a darker shade of gray than the shell; 

• each cell is composed - in this example - of 3 2 points. 

The larger dots on Fig. 1 represent the centers of the cells of size £ 0 = T 0 /3 = 3 admissible 
in the geometrical constructions referring to the cubes of such size. In this case, the minimal 
spacing between the centers of admissible cores equals 3. Considering L 0 = 9, we have the 
spacing £ 0 = L 0 /3. The admissible cells of a given size form a partition of Z d , and we denote 
by Cfc(x) the unique admissible cell of size £ k = L k /3, containing a given point x; naturally, 
Cfc(x) = C k (y) for all x,y from the same admissible cell, so there is no conflict with the 
previously introduced notation C k (c) where c was the cell’s center. 

It will be convenient to endow the set of the admissible cell centers c G C k in B Lk+1 (u ) with 
the natural graph structure, with edges formed by the pairs of nearest neighbors c, d with respect 
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to the max-distance, i.e., those with \c — c'\ = -4 = L k /3. Such a graph Bk+i will be called the 
skeleton graph of B Lfc+1 . Skeleton graphs will be used in Appendix A. 

The main tool for the analysis of the Green functions in such cubes often is the Geometric 
Resolvent Inequality (GRI). In its basic form, used in [18, 33, 19] and in numerous subsequent 
works, a single application of the GRI moves one from the center of a given cube B L (x) to (any) 
point y of the exterior boundary d + B L (x ) := {z : d(z, B L (x)) =1)}- Here d(-, •) stands for 
the graph-distance in the lattice Z d , with edges formed by the nearest neighbors in the norm | • |i. 
The notion of the exterior boundary is relative to an ambient set A D B L (x) (a subgraph of Z d ), 
when the analysis is carried out in a proper subset A of the lattice. As was said in Section 1.3, 
we employ another version of the GRI typical for the applications to the continuous Anderson 
models (cf., e.g., [26]). 

Given a finite subset A C Z d , we introduce the local Hamiltonian H\ := 1\H1\ [ f? 2 (A), 
acting in the finite-dimensional Hilbert space £ 2 (A) canonically injected into £ 2 (Z d ). II \ is self- 
adjoint; it is often considered as the restriction of H to the subset A with Dirichlet boundary 
conditions outside A, but the terminology here varies from one source to another (cf., e.g., [27, 
Sect. 5.2]). We introduce the interior boundary d~ A and exterior boundary <9 + A of the set A by 

<TA := {x £ A : dist(x, Z d \ A) = 1} , <9+A := d~ (Z d \ A) . 

For brevity, we often use notations like E A for the spectrum of II \, i.e., the set of its eigenvalues, 
counting multiplicity. In the case where A = Bl(u), we also write E UjR . This definition will be 
recalled where necessary. 

For the derivation of the GRI, it is convenient to use the language of the balls A r . of explicitly 
specified radius r, rather than the balls B R , with L directly related to the cardinality of each 
coordinate projection. Let A = A R k {c x ) C A', for some finite A' C Z d with d(<9+A, d~M) > 2, 
0 = l\ Rk l , and $ be the operator of multiplication by 0. Note that for any u £ A Rk l (c x ), i.e., 
any l u £ Ran <I>, we have the identities 1 A '1 U = l u = 1a1u» an d 

supp (Hl u ) C A Rk = A , 

since H is a finite-difference operator of order 2. Therefore, Hl u = 1 A IIl n . and similarly, 
using Ia'Ia = 1a, 

I^aI u 1a-^1a1u k R Hl u , 
so for any basis vector l u £ Ran <I>, 

Ha’^-u Ha\ u . 
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(2.3) 


As a result, one has the operator identity II\ <I> = H A i<&, thus for any E G E, 

(H a - E)$ = (H A i - E)$ 

= $(H a ,-E)-[$,(H a ,-E)}. 

Below, the energy E will be fixed and omitted from notation in the resolvents G A > = (7/y — 
E )~ l , G a = (H a — E ) _1 . Denoting 


W = [$,(H a ,-E)} = [$, A A z] (2.4) 

(note that [<f>, V — E] = 0) and multiplying the identity stemming from (2.3), 

$(H ai -E) = (H a - E)$ + W , 
by G a on the left and by G A f on the right, we obtain the identity 

+ G A WG A :. 

The RHS representation of W in (2.4) implies 


\\W\\ < 2||$|| ||Aa'|| < 8d. (2.5) 

Due to the relations (2.1) and (2.2) we have B Lk (c x ) = A R k (c x ), with R k = (. L k — l)/2. Let T k x 
be the boundary annulus of width 2 of B L k {c x ) (cf. Fig. 2): 

r‘ = lt := A fll (c r )\A Rl _ 2 ( C J, (2.6) 

and r^i = T k Cx (boldface notation) the operator of multiplication by l T k. More generally, in the 
case where R — (L — 1)/2, y G 7L d , B = B L (y), denote 

F| := A R (y)\A R - 2 (y), (2.7) 

and let r| be the operator of multiplication by l r |. Introduce also a shortcut for the minimal 
set of centers of L/,-admissiblc cells covering a given subset A C Z J : 

C k (A) = {ce C k : C fc (c) n A ^ 0}. (2.8) 

Observe that the Laplacian A, being the canonical graph Laplacian (on Z d ), annulates the con¬ 
stant functions, hence 

supp (A0) c r k Cx 

(in fact, supp (A0) = r k J. Indeed, for any 2 e A Rk (x) \ T k x , the function <ft = lA flfc _i takes 
the constant value 1 on the 1 -neighborhood of z, thus (A <p)(z) = 0. Similarly, 0 vanishes 
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supp 0 



Figure 2: Example for the GRI (2.12). Black dots on Fig.2 (a) represent the boundary of supp 0, and the white 
dots - the boundary of Z d \ supp 0. Their union is the set 10 = T* , shown as the dark gray annulus on Fig. 2 (b). 
The commutator W annulates any function supported by the complement Z d \ F0. 


on the 1-neighborhood of any point z e Z d \ A R k (c x ), hence (A 4>)(z) = 0. In other words, 
supp (A0) is covered by the union of the boundary of supp 0 = 1 Ah (c x ) and of the boundary 
of supp (1 — 0). 

It follows that the commutator W = [<&, f//y — E )] = [<f>, (A a > — E)\ satisfies the operator 
identity W = T k x WT k x , so for any subset AcA'\ A Rk (u), one has 

IaGvxI = 1 i'I'Cav!; + 1 aGvWGmI, 

= (UG^W(T k ,G kX i), 

and we come to the following form of the Geometric Resolvent Inequality: 

\\1aGa>xI\\ < \\W\\ • ||iAG A ,r*|| • \\r k x G AX i\\- ( 2 - 10 ) 


Introduce a slightly abusive but convenient notation, recalling that we are going to use a sequence 
of length scales following the recursion L = YkLk-i : 


11^ 


®L k (x)\ 


— C Wk \\l r kG 


B Lk ( x 


)(E) X k c 


a 


W,k 


_ \ r cL 

■— 1 k 


( 2 . 11 ) 


Here A symbolizes the decay from the center to the boundary of a cube. A more accurate (but 
cumbersome) notation would include the dependence of the symbol A upon the cube B. For 
brevity, let B' = B Lk+1 (u), B = B Lk (x). With A = r£ +1 (the set A appears in (2.9)-(2.10)), we 
infer from (2.10) and ||W|| = Y~ d Cwk 


\K +1 GB'xi\\<Y- d C w ,k • 

ll-rtfc+l/^ r\k 1 

II 1 u 

1 ■ ||r‘G B xt|| 

< Y d Cw,k • 

H r ‘G B x‘,|| 

E H r S +1 G , B'xS|| 


c: 

\C-Cx\=ik 

<||Gb|| A max Jlt +1 G B .4J|, 


( 2 . 12 ) 


12 












Figure 3: Example for the inequality (2.16). Black dots are the centers c £ C k (l'pJ of admissible cells covering 
the set Tg, shown as a dark gray annulus. Long thick arrows represent the terms r„ 1 Gb , Xc with c £ C /,: (l’g). In 
(2.17), x = u, so B and B' are concentric. 


as card{c : (c, c x ) G C k } < Y d . This bound is useful when ||Gb||^ is small. 
By self-adjointness of our Hamiltonians, we also have 


||Gb|| — Cw.fcUrBG'BXc^ll < Cw,fc||rB|| ||XcJ| (dist(E, E(iF B )) 

<Cw,k(dist(E,Y(H B ))~\ 

yielding an a priori bound useful in the case where 11 Gr 11 A is not small: 

Cw.k 


-l 


\K +1 Gb'Xc II < --rv max 

|u Cxl1 “ dist(£,£(f/ B )) MeC* 


-t/c+1 


Gb'XI 


More generally, in the case where (cf. Fig. 3) 

x G B = B L (w) C B Lfc+1 (u), dist(B L (u;),r^ +1 ) > 0, 


(2.13) 


(2.14) 


(2.15) 


with an arbitrary L < L^+i compatible with (2.15), the boundary belt T B figuring in (2.13) is 
covered by at most (in fact, less than) Y d cells C/-(c), c G C k . Therefore, in such a general case, 
we still have 


II < 


a 


W,k 


max 


ife+i 


dist (E,Y(H b )) cec fc (r|) 


Gb'Xc 


(2.16) 


In Appendix A, we will use (2.16) in the situation where x = u, thus B = Bl(u) is concen¬ 
tric with B' = B L k+1 (u), and dist(£’, E(iT B )) = : £b > 0, so 


-t/c+1 


G 


Br (u)Xr II < -- maX 

B Lk+1 (u)A Cx \\ - £b cgCfc(r , } 


a 


W.k 


-ifc+1 


Gb'Xc 


(2.17) 
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3 Dominated decay and EVC bounds. Main results 


Below the centers of various cubes B,(«) will be assumed to be the admissible centers at the 
respective scale. 

Consider a cube B = B Lk+1 (u) along with its skeleton graph B defined in Section 2 (with 
the vertex set B fl C k ), and introduce the function / : B —* M + given by 

f : x ^ \\r k u +1 G BX k x \\ . 


Then by GRI(2.12), 

f( x ) < ||G Bifc( *)|| A max /(c) . ( 3 . 1 ) 

\C)X f tc 

An inequality of the form (3.1) is most useful when 11 (a;) 11 ^ < q < 1; in this case, using 

an iterated application of the GRI, it is not difficult to prove the bound f(u) < q Y ~ 1 . Below we 
formulate a more advanced analog of this simple bound (Lemma 1) in a more general situation 
where for some S > 1, there are at most S vertices c G B where ||B ifc (c)|| A fails to be small. 
Recall that stands for the spectrum of H-q l ^ u ), and Y k is an odd integer, so 1/ = 2A/ + 1. 

Definition 1. Let be given an integer k > 0 and real numbers e > 0 and E. 

• A cube B L (u) is called ( E , e)-NR (non-resonant), iff dist (£ Uji , E) > e; 

• A cube B Lk+1 (u) is ccdled (E, e)-CNR (completely non-resonant), iff for all j — K k ,, K k + 
Yfc+i — 2 the cube BjL k / 3 (u) is ( E , e)-NR. 

The role of the concentric cubes B jL k / 3 (u) will become clear in Appendix A, in the proof of 
Lemma 1 (cf. also Fig. 4). 

Definition 2. Let be given an integer k > 0 and reed numbers e > 0 and E. A cube B L k (v) 
is ccdled ( E,e)-NS (non-singular), if E f E , i;L , and |G'h,. a ( w ) || A < e. Otherwise, it is called 
(E, e)-S (singular). 

Below we choose the sizes L of cubes Bffu) and the parameter £ > 0 figuring in Definitions 
1 and 2 in a specific way. First, we take L e { L k , k > 0}, with L k = Y k L k _i and Y k defined 
in (4.7); £ = L^ bk in the context of Definition 1, while in the property (E,e )-CNR we set 
£ = L~)( Sk , with recursively constructed sequences ( b k , k > 0) and ( s k , k > 0) (cf. Table (4.14)). 
Now we are ready to formulate the main results of the paper. 


Theorem 3.1. Consider the random Hamiltonian H(u) of the form (1.1) and assume that the 
marginal PDF F v of an IID random potential V( -: ut) is Holder-continuous of some order (3 G 
(0,1]. Further, assume that for some b 0 > d//3 and L 0 G N satisfying 


L 0 > L 0 (rj,T ) 


max 


11 1/t2 , 9 


4(6d+77) 



8 " 



(3.2) 
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Figure 4: A 1-dimensional example for Definition 1. Here I\ k = 3, Yfc+i = 2 ■ K k + 1 = 7. The thick black 
interval shows the position of the core of the cube Bi fc+1 (0), with some k > 0 and L k +i = Yk+iLk, decomposed 
into cells of size £ k = Lfc/3 (short gray intervals). The CNR property requires the Yk+i — 1 = 6 long black 
intervals (cubes, for d > 1), to be non-resonant. These cubes may be used in the geometrical procedure in the proof 
of Lemma 1 (cf. Appendix A), along with the cells shown as thick gray intervals. Thinner gray intervals are the 
remaining (^-cells. 


with t = and rj = |(/ 3b 0 — d) (cf. (4.1)), the following condition is fulfilled: 

P {B Lo ( 0) is (E, Lf bo )- S} < (3Y] - 4)" 2d = 529~ d . 

Then there exist positive sequences (8k)k>i, {Kk)k>i such that lim 6k— lim n k = land for 

fc—>+oo k —^-)-oo 

all k > 0, 

sup p{b L;c is <e~( Lk)Kk . 

See Section 5 for explicit bounds on 8 k and K k . 

Let be the set of all bounded Borel functions f : M —» C with ||0||oo < 1. By well-known 
techniques (cf. Sect. 6), one can infer from Theorem 3.1 the following result. 


Theorem 3.2. Under the hypotheses of Theorem 3.1, there is a function f : R + -+ R + which 
admits exponential scaling limit, viz. 


lim 

L —l~oo 


lii In/(L) 
In L 


1 , 


(3.3) 


and such that for all x. y + 7L d and for any sufficiently large domain AC 7l d containing x and 
y, one has 


E 


sup 

4>e&i 





< f([x-y |) • 


(3.4) 


The principal analytic tool used in the proof of Theorem 3.1 is the following 


Lemma 1. Let the integer sequences ( Y k ), {S k ), (L k ) and positive real sequences ( b k ), (s k ) be 
defined as in (4.13)-(4.14). Fix some k > 0 and suppose that a cube B Lk+1 (u) 
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• is (E, L k Sk )-CNR, and 

• contains no collection of(S k +i + 1) pairwise disjoint (E, L(( l>k )-S cubes of size L k with 


admissible centers c G C k . 


Then one has 


\G BLwM ( E) r<LfL 


d/s T -b k (Y k+ i— 55fc + i—1) 
'k L k 


(3-5) 


See the proof in Appendix A; like its counterpart from [26], it is in essence a variant (or 
rather an adaptation) of a well-known argument going back to [19]. 

As usual in the MSA, we also need an eigenvalue concentration (EVC) estimate to bound 
the norm of the resolvent near the spectrum. 

Lemma 2. Assume that the marginal probability distribution of an IID random potential V is 
Holder-continuous of order [3 G (0,1]. Then for any cube of size L one has 


P {B L (u) is not (E, L“ S )-NR} < Const L d L ~ ps . 


(3.6) 


In the case where V admits a bounded probability density, hence (3=1, this is the classical 
result by Wegner [34]; see also a short proof, e.g., in [15]. A simple adaptation to Holder- 
continuous (and more generally, continuous) marginal distributions, sufficient for our purposes, 
can be found in [12], where it is shown that an EVC bound for the potentials with Lipschitz- 
continuous marginal PDF F v can be automatically transformed into its counterpart for the PDF 
with an arbitrary continuity modulus. Optimal Wegner bounds have been proved earlier for 
various types of operators; cf., e.g., [9, 10, 29]. 

4 Adaptive feedback scaling 

4.1 Technical assumptions and some useful inequalities 

In the recursive construction of the sequences (b k ) k > o and (s k ) k > 0 , mentioned in the previous 
section, the crucial input parameter is b 0 . Given the marginal distribution F v of the random 
potential V : 7L d x f) —* R, which we assume Holder-continuous of order (3 G (0,1] until 
Section 7, we always assume that b 0 > d/(3 and introduce the scaling parameters 


V ■= ^{/3b 0 - d) > 0, 


(4.1) 



(4.2) 
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The initial length scale L 0 is always assumed to be large enough, to satisfy the explicit condition 
(3.2). Further, set 

*i = 9, Si = 1, a 1 = (3y 1 -4) d . (4.3) 

We make a crucial assumption, 

Po ■= P { B Lo ( 0) is not (E, L~ b °)- NS} < a/ 2d = (4.4) 

and introduce the parameters 0 O e (0, 1/3) and a 0 > 0 by letting 


lnai l + 30o In po * 

In pQ 1 2 ’ 0 In L 0 


(4.5) 


The scale-free probability threshold in the RHS of (4.4) is slightly better than 841 -d given in 
[26] . This marginal modification is due to a geometrical strategy of the proof of Lemma 1 which 
deviates from that of an analogous argument in [26]. It is clear, however, that the importance of 
the scale-free probability bounds from [26] goes far beyond the explicit numerical estimates for 
specific lattices. 

Further, introduce an integer 


ft = ft(po, *i) := min {k > 1 : (1 + 6 0 ) k > 2 d/a 0 }, (4.6) 


and define the integer sequences (Y k ) k > i, (S k ) k >i, and (L k ) k >\ as follows: 


Y, = 


S k := 


Next, let 

L k ~\-\ ■ 1 fc- 

and define for k > 0 (cf. (4.2)) 


n = 9, 

k< ft, 

T 1/8 
L k- 1 

5 

k > ft, 

fSi = l, 

k < ft, 


k > ft. 

k +1 -^fc+l I? 


4 1 

b k + i = ~N k b k , Sfc+i = 2^ fc+1 ' 


(4.7) 

(4.8) 


(4.9) 


(4.10) 


4.2 Main formulae and relations between technical parameters 

The probability p 0 is small enough, viz. 

p 0 :=P{B i0 (0)isnot(T;,L- fe0 )-NS} <23" m , (4.11) 
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and the initial length scale L 0 is large enough; specifically, it suffices that 


Lq > L 0 (rj,T ) := max 


1 / .2 4(6d+7j) — — #■ 

ll 1/r , 9 , p 0 3 ^, p 0 b ° 


(4.12) 


It is to be emphasized that the parameters r, rj, p 0 figuring in the hypothesis (4.12) are related 
only to d and b 0 , so there is no vicious circle in (4.12). 


1 = Si < Yi = 9 

Ni := Y x - 5S 1 - 1 = 3 

A 

o 

-O 

77 = |(/36 0 — d) = f3s 0 — d > 0 

ai = (3Yi - 4) d = 23 d 

0 i= H T = ik 

1 lnai 1+36 »q a ^ 1 

1 - 2 > 4 < 3 

•8 = min{fc: (1 + 6 0 ) k > jg} 

(To := hi /Iq 1 / In L v 

- — InVi 30 o „ 

T ° mm |_lnL 0 ’ l+30o ’ T \ 


j Y\ — 9, k < R 

k+1 ~ \l L l\ > 

/Si = l, 

+1 Uin+ij, *>* 

N k = Y k -5S k -l>3 

A* = (|) fc 7V 1 ---7V fc 

b k = A k b 0 = | s k , k > 1 

a fc = (3Yfc - 4) d 

^ _ f 1 + 4 ^ (l) |) ; k <&+ 1 

|(iS , fe + 1), k > A + 1 

&k = B fc a 0 — B k - ■ ■ B 1 cr 0 

\|(5 fc + l), k>& +1 

£fc= ■ £>2 £> 1 , k > 2 

4 > ln(5fc + 1) ” ln(3/2) / 1 

Pfc < Or k <Q k 

P*:=P{B t .(«) is (B,L^)-S} 

roA- < 

tDfc := IP {B Lk (u) is not (E, L^)-CNR } 


4.3 Unbounded growth of the geometric scaling parameters 

The following statement is an important ingredient of the proof of exponential scaling limit in 
the scheme with varying scaling parameters Y k , S k (cf. Sect. 5). 
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Lemma 3. Let be given an integer L () > Lo(rj), with L (] (//) given by (3.2). Then the sequences 
( Lk ) /,:>o, {Y k ) k> si and (S k ) k> & given by (4.7)-(4.8), are strictly monotone increasing, and for 
all k > 1 one has 

^Y k <S k < ^Y k . (4.15) 

Proof. Let L 0 > 11 1</t2 , k > &. Then we have r k = r and 


n+i = \li\ = [YfLui = L(L lu\) t lu\ > L(L^SJ) T ^*-iJ 


> 

> 


( ir /r2 ) 

io r/T ' 


J k -1 


> 


( 10 t/t2 + 1 ) 


J k-l 


k -1 


> 10 L Ll_,\ = 10 Y k > Y k . 


Furthermore, 


S k+ i — 


1 



1 


1 

-Y k+ 1 
9 + 

> 

10 

-Yk 

9 

> 10 

-Y k 

9 


> s k . 


Therefore, the sequences (Y k ) k> % and (. S k ) k> & are strictly increasing. 

To prove the LHS inequality in (4.15), notice that Y &+1 > 10Y$ >10-9, and for any real 
2 / > 90 one has [|J > | — 1 > -jfc, hence for all k > & + 1 


10 


Y k <S k = 


-Y k 


9 


< -Y k , 
~ 9 ’ 


as asserted. 


(4.16) 

□ 


4.4 Scaling of the GFs 


Lemma 4 stated below is the usual analytic component of the MSA, providing the conditions 
under which localization bounds at a scale L k imply similar (or better, in this case) bounds at 
the next scale L k+1 . Its proof is very short, but the bulk of technical work is actually done in b 
Lemma 1 proved in Appendix A. 

Lemma 4. Define the integer sequences (Y)), (Sf), ( Lj ), (Nf and positive real sequences (b :j ), 
(. Sj ) as in (4. 13)— (4. 14). Assume that a cube B Lk+1 (u), k > 0, is (E , Ljjjf)-CNR and contains 
no collection of S k+] + 1 pairwise disjoint (E, Lf bk )-S cubes of radius L k with admissible 
centers. Then the cube B L k+ 1 (u) is (E, L~ bk+1 )- NS. 

Proof. By Lemma 1, we have ||G Bi (n)|| A < L - bk N k+1+ i 

Therefore, recalling that b k > 

d/(3 > d > 1, 24+1 > 3, we obtain 


In G 


BL fc+1 («)| 


In L k 


> bkNk+i ~ — > b k N k+ 1 



23 

> ^ k 2A^ k+1 
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With L k+ 1 = L k Y k+ 1 , In Y k+l / In L k < 1/(16 d) < 1/16, we have 


b k -^jN k+ 1 • -— > ~N k+ i b k — b k +i ■ 
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□ 


4.5 Scaling of the probabilities 

In the next statement, we establish an important technical ingredient of the proof of the key 
Lemma 6: an upper bound on the probability of “tolerated resonances” following from the 
Wegner estimate. As usual in the MSA, such upper bounds essentially shape those on probability 


of “insufficient” decay of the Green functions and, ultimately, of the eigenfunction correlators. 


Speaking informally, one cannot get bounds better than those stemming from a Wegner-type 
analysis, so we have to make sure the latter is compatible with the exponential scaling limit. 

Lemma 5. Define the sequences of positive integers ( L k ) k > 0 , {Y k ) k > 1; ( S k ) k >i and positive real 

12cH-4 

sequances (tt> k ), ( q k ), (g k ) as in (4.13)-(4.14). Assume that Lq > 9 ^ (cfi (4.12)) Then the 
following bound holds 3 ; 


Vk> 1 W k <q k :=2L k Bk . 


(4.17) 


Proof The probability ro ^+1 refers to the event (cf. Definition 1) that for at least one integer 
j G [K k , K k +Y k+ i — 2], the cube B jLfc / 3 (w) is (E, e)-R. The probability of each of these Yfc+i —1 
events is bounded by Wegner estimate (cf. (3.6)), thus 



(4.18) 


By hypothesis, L 0 > 9 n > e /v , so for k = 0, (4.18) becomes 


12d+4ri 



(4.19) 


the last RHS equality is the definition of Qi in (4.13). Hence < L 1 ei = \q x . 
Now let k > 1. By (4.14), (3s k > fisi = = 2 /3b 0 > 2d > 2, hence 


_ 5 4JVi 



Using the factor ”2” in the RHS of (4.17) might seem artificial, but it becomes convenient in (4.26). 
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where 


1 In Y k+1 1 In Ll /{Wd) 1 
A,. >i- fe+1 >_ h _> _ 

' 2 2 In Lfc+i “2 21nL fc 4’ 

so rcife+i < 2L /: = 2L^j +1 , where (oj) follows the recursion 


Qk+l 

Qk 


i>k 


bk-i 


—— > ~r(Sk + 1) > — (<Sfc + 1) — -Dfc+i 

5 5 6 


(4.20) 


Since toi < qi — L l 81 , this implies t Vj < qj = L ■ Sj for all j > 1, owing to the definition of 
Qj= Dj■■■ D 2 Ql (cf. Table (4.14)). □ 

Recall that we defined in (4.6) an integer if = min{/c > 1 : (1 + 0 o ) k > 22}. Now we turn 
to the heart of the adaptive feedback scaling analysis. 

Lemma 6. Consider the sequences of positive integers ( L k ) k > 0 , (Y k ) k > 1; (S k ) k > i defined as in 

12d+4?7 

(4.7)-(4.8). Let { p k , k > 0} be defined as in LHS equation (4.25). Assume that L 0 > 9 ^ 
(cf. (4.12)), and one has (cf. (4.11)) 

p 0 < (3Yi - 4)~ 2d = 23 _m . (4.21) 


Define recursively a sequence of positive numbers (a k ) k >o ■ 

(To = &o = lnp^/lnLo, 

(Tk = B fc <7 0 , B k = Bi - ■ ■ B k: 

jj _ f 1 + $0 £ (l, f)j j = l,...& 

\§(*Sj + 1) — j > A.+ 1. 

Then for all k > 1, the following bound holds: 

Pk = p{b L i («) is S } =: Lf- < Lp- . 


(4.22) 

(4.23) 

(4.24) 


(4.25) 


Proof. By Lemma 4, if B Lk+1 (u) is (E, L~ k + k f 1 )-S, then either it contains at least (S k +i + 1) 
pairwise disjoint ( E , L ^ bk )-S cubes of size L k , or it is not ( E , L^)-CNR. Since S k + 1 +1 > 2, 
it follows by a simple combinatorial calculation that, with a k+ 1 = (3Y k+ i — 4) d (cf. (4.14)), 


Pk+i < 7^{ a k+iPk) Sk+1+1 + (4-26) 

By Lemma 5, we have tv k+ i < q k+ 1 = L k +\ +1 , thus 

Pk+i < a k+1 p k ) Sk+1+1 + ^q k+ 1. (4.27) 
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Now we shall analyze separately several cases: (Al), (A2), (Bl), and (B2). 

(Al) Suppose that for some k Q 6 [0,.fi — 1] and all k e [0, k Q ] one has p k +i > qk+i • ( The 
analysis of the case (A2) below shows that once Pk+i < qk+i, the inequality pj < q,j then 
follows by induction for all j 6 [k + 1, ft].) Observe that 

^(a k+1 pk) Sk+1+1 > Pk+i ~ i qk+i > ~Pk+1, (4.28) 

so pk+i < (ak+ip k ) Sk+1+1 ■ By finite induction in A; e {0,..., k Q } (which we are carrying out 
now), we know that d k = In p^ 1 / In L k > cr 0 : this is true for k — 0 due to (4.22), and for larger 
k this will stem from the inequalities (4.30)-(4.33) we are going to prove. Therefore, starting 
with k — 0, we can write 


^ dlnafe + i^ dlna k+ i > din ax din a! 1 + 39 0 

dfclnLfc d k In L k ~ a 0 In L 0 lnpo 1 2 

and with S k + i + 1 > 2, In L k+ 1 = In L k + In Y k + 1 , Y k + 1 = Y\ , we obtain 


In Pfc+i 

In L k+ 1 


&k+l L &k ' (‘S'fc+l + 1) ( 1 — 


d In dfc+i \ In L k 


>a k -( 1 + 3d 0 ) • 1 


efic In L k y In Lfc+i 
In Y, 


In L k + In Yi / 

> Ok ■ (1 + 3&o) ■ (1 - To) (see (★) below) 

> cr fc • (1 + d 0 ) = = (7 fc+ i (cf. (4.24)), 


(4.29) 


(4.30) 

(4.31) 

(4.32) 

(4.33) 


so as long as p k+ 1 > q k + 1 , we obtain inductively that p fc+1 < L fc _^ +1 , with 

4 

o'fc+i = -B/c+iOT, -Bfc+i = 1 + 6*o < - = D k +\ ■ 

(★) To derive (4.33) from (4.32) we used the condition r 0 < from (4.13). 
Thus (4.29), (4.30)-(4.33) hold inductively for all 0 < k < k D , as claimed. 

For further use, we need to make the following observations. 

♦ By (4.5), a 0 := and by (4.13), q x = thus 


ai = B x oo = (1 + 6* 0 )o- 0 < -(Tq 


4 Infa) 1 ) 

3 In L 0 



_ 8 _ 

where the last inequality is equivalent to the assumption fro > Po 3v (cf. (4.12)). 
♦ We also have cr 2 <Q 2 , since 


o o 16 

C2 — B^BiOq < —Co, 

y 



12 

8Y3 


b 0 = 
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so a 2 <Q 2 stems from 


which, in turn, follows from the assumption L 0 > p 0 b °, equivalent to lnp 0 1 / In L 0 = a 0 < | b 0 
(cf. (4.13)). 

♦ More generally, the inequality 

cTj < Qj (4.34) 

follows inductively from cr 2 < g 2 by the recursion a 3 = B 3 a 3 - i, g 3 = I) 3 o 3 -i, as long as 
Bj < D 3 , and the latter holds by (4.14): 

Bj = 1 + 6 q < - = Dj, j < if, 

Bj = - (S 3 + 1) = D 3 , j > .ft . 

The conclusion of this step of analysis is that for all k e [0, k 0 \, 

Pk < Ll ak , a 0 <cr 1 <a k < g k ■ 


(A2) Consider the case where, for some k < .ft, 


j-(T k+1 _ , _ j-Bk+1 

Ll k +1 — Pk+l Y Qk+l ~ L k+l 

First, let us show that in this case (a k+2 p k+ 1 ) 2 < q k +2 , so by recursion (4.27) one has p 3 < q 3 
for all j = k + 1,..., .ft. Indeed, with Y' k+] = V) = 3 2 , 1 = 1 and 

ai = (3Yi - 4) d < 3% d = Y; id/2 , L k+2 = YiL k+ i, 

the required inequality would stem from (Y^ 2 p k+ ij < (YiL k+ i)~ 0k+2 . We have assumed 
Pk+i < q k + i = L^ Qk+1 , so it would suffice that 

\y 2>d T ~‘2‘Qk +1 ^ \r — Qk+2 T ~Qk+2 

J 1 ^k+l — r l ^k+l J 


be., (2 g k+ i — Q k+2 ) lnL k+ i > (g k +i +3 d) In Y x . Since g k + 2 — \ g k +i, hence 2 g k+ 1 — g k +2— 

- g k+ 1 , it suffices to validate the inequality 


equivalent to 


g k + 1 hi L k+ \ > (g k+ 1 +3 d) In Y, , 


1 nL t+1 > 4(gw + M > i nFl =4 bl + — I W, 


gk +1 


^fc+i 


(4.35) 
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4(6d+rj) 


Since L 3 > L 0 , pj > p k = rj/2, the relations (4.35) stem from the hypothesis L 0 > Y 1 v = 

4(6d+r l) 

9 'i , by a straightforward calculation. As was said, by (4.27) this implies 

V j G (1,..., L~ aj = p 3 < q 3 = L~ ej < LJ aj , as a 3 < o 3 . (4.36) 

The logic of the analysis for k > & is quite similar to that of the cases (A1)-(A2), but here 
one has to operate with growing sequences (Y k ), (S k ). 

(Bl) Suppose that for some k Q > R and all k E [£, k Q ], one has p k+ 1 > q k +i- By the same 
argument as in (Al) (cf. (4.28)), it follows that p k+k < (a k+ ip k ) Sk+1+1 . Unlike the case (Al), 
now = B k - ■ ■ Bi > 2d/a 0 , by definition of hence 


a k = Ba-Uq > 2d. 


As a preparation for the bounds (4.37)-(4.41), recall Y k+1 = [L T k k+1 J < L T k +1 , r k+ 1 = l/(16d), 
3Y k+ i < Y k+1 , so it follows that 


In ((3Yfc+i) d ) 2d In Y k+1 < 2dr k+ i In L k < 1 
crfcln L k a k In L k ~ a 0 B k In L k ~ 8 

since coBa > cr 0 B.i ? > cr 0 (l + > 2 d (cf. (4.6)). Therefore, 

InPfe+r ^ , 0 , lN ln(a fe+ ip fc ) lnL fc 


In L 


k -\-1 


— O'k+l > (S k+1 + 1)- 


Y a k - {S k+ 1 + 1) ( 1 


In L k In L k+ i 
2d In Y k+ 1 


a k In L k J In L 


In L k 


jfc+i 


> (Tfc • {S k +1 + 1) • ( 1 — — 


In Li 


8 / In L k (l + T fc+ i) 

^ , 7 16 

> a k ■ (S k+ 1 + - 1 -) ' g ‘ Y 7 

2 

> a k ■ -(*S'fc +i + 1) > B k +l cT k ■ 

Further, N 3 = Y 3 — 5 S 3 — 1 > 9 S 3 — 6,5' ; = 3 Sj, and 1 + 0 (] <4/3, thus 


AM. 12 14 

D j = —- > —g > —Sj > max 
5 5 9 


2 (Sj + 1 ), 1 + 6*0 


— Bj. 


(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 


Remark 1. A subtle (albeit merely technical) point here is the definition of B 3 and D 3 at j = 
^ + 1, which may have attracted the reader’s attention in Table (4.14) and in (4.24): switching 
from the initial values, resp., 1 + 9 0 and 4/3, to the large ones growing with j, is delayed by 
one step and occurs at j = ^ + 1. This is due to a one-step delay in the multiplicative recursion 
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Qj + 1 = Dj + 1 Qj= -y- Qj in Lemma 5 (cf. (4.20)). Here we rely on the value Nj, not Nj + As 
such, this fact is not directly related to the relation cr$ + i = B& +1 a&, but it is convenient to keep 
the general bound g 3 > a 3 , and this is why we set, artificially, B & + 1 = B& = 1 + 8 0 , so as to 
have Am = | > Am- 

Now observe that Qk+i> <?k+i for all k, since we have (cf. (4.20)) 

Qk+i _ D k Q k > Qk > y Qi y ^ 
o'/c+i B k (Jk <Jk (J\ 

(B2) Suppose that for some k > R, one has p k +i < Qk+ i- To show by induction that pj < q 3 for 
all j > k + 1, we need to check, starting with j — k + 1, that 

< Qj+i, (4-42) 


for (4.42) would yield immediately 

Pj +1 < \ ( a jQj) Sj+1 + \qs +i < \oj +1 + 7^+1 = 9 J+ i = Am +1 . 

The relation (4.42) is equivalent to 

(3 2d Yf d y 


(4.43) 


j2d\/-2d\ (Sj+l)e: 


'j v W + 1 r W + 1 

- Am A 


Further, on account of 3 Yj < Yp /2 (as Y, > Y\ = 3 A it suffices that 


Y^d(Sj+l) (Sj+l)gj ^ y Bj + 1 L ^ + 1 

Recall Yi = L]^ for i > A and o J+ i = Dj + i Qj. By direct inspection, each of the inequalities 
(4.44)-(4.45) follows from the next one in (4.44)-(4-.46): 


jpSj+\)Qj —Qj+i ^ ^3d(Sj+l)+gj+i 

gj ^ ^ ft' 


(A + 1) ( T — 


16 


> 1 + 


16 


D 


3 +1 


(4.44) 

(4.45) 

(4.46) 


Finally the validity of (4.46) follows from D j+1 = | (S' J+ i +1)• This proves (4.42), hence (4.43). 


The assertion of the lemma is proved. 


□ 


♦ This marks the end of the ’’renormalization group” analysis of the Green functions. In 
the next Section 5 we interpret the finial outcome of this analysis as asymptotically exponential 
decay of the GFs, and then derive in Section 6 similar decay properties of the eigenfunction 
correlators. 
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5 Exponential scaling limit (ESL). Proof of Theorem 3.1 


We have proved that for all k > 0, 

<L~ k ’\ (5.1) 

The language of “polynomial” bounds - with growing exponents - has been so far convenient 
but certainly looked rather artificial, so our next goal is to show that (5.1) can be interpreted as 
follows: 

p{b l , is (£,e- (ifc) 5fc )-s} < e~ {Lk)K \ 

where S k , n k /*• 1 as k —* +oo. This is a matter of simple calculations. 

Indeed, by induction, b k = A k b 0 . Since Sj < Yj/ 9, we have 


Nj > Y j ~ 5Sj - 1 > 3Sj > -Yj. 


Therefore, 


4 > M-‘n y j =-r 3 "‘ L ‘. 


3 =1 


1—r- 7 — (in 1- 

_ £ lni fe V b o ) ^ T 


In Ln+2fc 


> L, 


k k 

since Yj /* +oo, thus k/ In L k —y 0. Consequently, 


= L 


i-o(i) 


7 , r rl-o(l) T l-o(l) 

~ b k << ^-lnL fc -L. _ „~c k L k 


L~“ k < e 


= e 


C k —> + oo. 

k ->-|-oo 


More precisely, L k ~ ^^+ 1 )/( 16d ) f or k ^ so f or some 1 < q ^ 1M±1 


16 d 


(5.2) 


(5.3) 


In L k >C + C'q k ~ A > C"q k . 


Thus 

In In L~ h C'" 

_ h _>i_ 

In Lk ~ (l + e) fc ’ 

Similarly, for the probabilities p k < L]~° k we have 


£ « 


1 

16 d 


In/ 4 1 > a k In L k > a 0 5i... B k , 


where Bj ^ CSj ^ C C 0, for 3.11 3 > Y. By taking a sufficiently small constant 
C" > 0, one can extend this lower bound to Si,..., B&. 

I1174 1 > C"c k Y { • • • Y k > C"'c k L k > Ll~ a(k \ 

with a(k) < h k , h G (0,1). 
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6 ESL for the eigenfunctions and their correlators 


It is well-known by now that a sufficiently fast decay of the Green functions, proved with suf¬ 
ficiently high probability at each energy E in a given interval /CM, implies both spectral 
localization (a.s. pure point spectrum in / with rapidly decaying eigenfunctions) and strong 
dynamical localization, with rapidly decaying averaged EF correlators. Such implications can 
be established with the help of different methods. For example, in the bootstrap method pre¬ 
sented in Ref. [26], the fixed-energy estimates in probability, proved at a given energy E 0 , are 
extended to an interval I 0 = [E 0 — e, E 0 + e) with sufficiently small e > 0, by means of the 
energy-interval (a.k.a. variable-energy) MSA induction; the core procedure goes back to earlier 
works [18, 33, 19]. 

In an earlier work [13] (cf. also the book [16]) we proposed an alternative approach based on 
an argument employed by Elgart et al. [21] in the general context of the FMM and encapsulated 
in a fairly general, abstract spectral reduction (FEMSA => VEMSA). Similar ideas, in essence 
going back to the work by Martinelli and Scoppola [30], were used in other papers; cf., e.g., [ 8 ]. 
Introduce the following notation: 

F XjL (E,uj)= max \G BL ( x )(x,y; E,uj)\. (6.1) 

yed~B L (x) 

We formulate the spectral reduction in the following way (cf. [13, 16]). ( Notice that the 
quantities b/, are unrelated to the sequence of scaling exponents hi,..) 

Theorem 6.1. Let be given a bounded interx’al I C M, an integer L > 0, two disjoint cubes 
B L (x), B l{u), an d the positive numbers a/,, b^, c l, q l satisfying 

b L < min [a L c 2 L , c L \ (6.2) 

and such that 

V E e I max P{F 2 , L > a L } < q L . 

ze{x,y} 

Assume also that, for some function f : (0,1] —$■ M + , 

Ve G (0,1] P {dist (£(// Bi( x)), X(H BL(y) )) <e}< f(s ). (6.3) 

Then 

P (sup max F Z . L (E) > a L | < + f{2c L ). (6.4) 

(Eci z £{ x ,y} J dl 

Consequently, taking into account the results of Section 4, for some (f - /* 1 as k —> +oo, one 
has 

P(3 Eel: B l (x) andB L (y) are {E,Lf bk )-S | < e~ L ' k . (6.5) 
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The proof given below is based on the following 
Lemma 7. Let be given L > 0 and positive numbers o l, b^, c l, q l such that 

b L < min \a L c 2 L , c L \ , (6.6) 

supP{F 2iL (£) > a L } < q L . (6.7) 

Ed 

There exists an event B z such that P {B z } < b^q^ and for any uj ^ B z , the set <§ z (f2 a) := {E : 
F z .l(E) > 2 a r .} is contained in a union of intervals 


M z 

U Ij, Ij '■= {E : \E — Ej\ < 2c l } , M z <(3L) d , 

3 =1 

centered at the eigenvalues Ej G P ( //(>')) ft I. 

Proof Consider the random subsets of the interval / parameterized by a' > 0, 

g{a'\u) = {E: F Z , L (E) > a'} 
and the events parameterized by b' > 0, 


B(b') = (cc G U : mes(<f(a) > b'} = < u G f) : / If z<L (E)>a L dE > b' > . 


Using the hypotheses (6.6)-(6.7), apply Chebyshev’s inequality and the Fubini theorem: 


P {B(b L )} < b^E[mes(^(a L ))] 

= b L 1 J dE E [ If z , L {E)>a L ] < t>^ 1 P {F z ,l(E) > a L } . 

Fix any u £ B(b L ), so mes(^(a L ; cu)) < b L . 

Further, consider the random sets parameterized by d > 0, 


H(c') = {A G M : min |Aj(cc) — A| > c'}. 

3 

Note thatfor a L G (0, c l ),A^ l := {E : dist(iP, 1Z(2c l )) < b L } C 1Z (c L ),hence A c bL = M\d bi 
is a union of intervals at distance at least c l from the spectrum. 

Let us show by contraposition that, for any u> £ B(bL ), one has 


{E : F z ^(E]uj) > 2a/,} P\7Z(2cl) — 0. 
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Assume otherwise, and pick any point A* from the non-empty intersection in the LHS. Let 
J := {E' : \E — A*| < b L } C A bL C 1Z(C). By the first resolvent identity 

\\Gb l (z)(E')\\ > ||G BlW (A*)|| - I E' - A*| \\G Bl{z) (E')\\ ||G Bi;W (A*)|| 

> 2a l - b L • (2c l) _1 (cl) _1 > a L , 


where in the last line we used the assumption (6.2). We also used the bounds ||Gb l (z)(A*)|| < 
(2 Cl)- 1 and 


\\Gbu*)(E')\\ < 


< 


< 


dist (-£7', S) dist(A*, S) — \E' — A*| 2cl — be 


where b L < c L . Consequently, the entire interval (A* — b L , X* + b L ) of length 2 a L > b L is a 
subset of <§{a_ l;cu), which is impossible for any u jL Bib[j). This contradiction completes the 
proof. □ 

For ctL, Cl < 1, which is a frequent situation where Lemma 7 is applied, one can give a 
simpler variant of the bound (6.4). 


Corollary 1. Let be given a bounded inten’al I C M, an integer L > 0 and disjoint cubes 
B l (x), B l{u)- Assume that for some a^,, q l G (0,1] 


sup max P{F Z , L {E) > a L } < q L . 

Eel z&{x,y} 

Assume also that, for some function f : (0,1] —> M + , 


( 6 . 8 ) 


Ve e (0,1] P {dist (S (H Bl{x) ), UH PlL(y) )) <e}< f(e) (6.9) 


Then 


P < sup max F ZjL (E) > max 

l Eel ze{x,y} 


1/2 

ClL, q L 


< W q T + /(2qf) . 


Proof. Let a’ L = max [a L , q/ 2 ], c L = q/ 1 ,then tv. := a,cj = min[aic|,c t ], since at. c, < 1. 
Thus (a' L ,bL,CL) fulfill the condition (6.2). Further, the function s i-)- P{F, L (i7) > s} is 
monotone decreasing, and a' L > a L , so we have P{F, iL > a' L } < t\ L , and the claim follows 
from (6.4). □ 


Proof of Theorem 6.1 Define the events B x . B y related to the points x, y in the same way as 
the event B z relative to z in the proof of Lemma 7, and let B = B x U B y . Let x f B. Then 
for both values of z E {x, y}, the set <f z (a) is contained in the union of the intervals J Z:i = 
[E- 2cl, e\ z) + 2c L ]. Therefore, 


P 


c a: inf max F - AE) 
Eeize{x,y} ’ V 


> a L 


P{(U : dist(S x , 'Ey) < 4c L ); 
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the latter probability is bounded with the help of (6.3). 


□ 


Now the derivation of strong dynamical localization from the VEMSA-type estimates (6.5) 
can be obtained in the same way as in [26], directly in the entire lattice Z' :/ . This requires an a 
priori polynomial bound of Shnol-Simon type (cf., e.g., [31, 32]) on the growth rate of spectrally 
a.e. generalized eigenfunction; the latter becomes unnecessary in arbitrarily large finite balls (cf. 
[13, 16, 17]). 

Theorem 6.2 (Cf. [13, Theorem 7] ). Assume that the following bound holds true for some 
£ > 0, L G N and a pair of disjoint cubes B L (x), B L (y): 


P{3 E G I : B l {x) and B L (y) are (77, e)-S } < 


Then for any cube B L fw) D (B i+1 (x) U B i+1 (x)) one has 


E [ | (lx | 0 (77 b) | lj,|] < 4e + f Jl- 


( 6 . 10 ) 


A 


We had incorporated in the definition of the non-singular cubes, through the ’’norm” • , 

a combinatorial factor measuring the volume of a cube 11 / ( •); otherwise, such a volume factor 
would be present in the term 4e in the RHS of (6.10). 

The extension of the EFC decay bounds to the entire lattice can be done with the help of the 
Fatou lemma on convergent measures; such a path was laid down in earlier works by Aizenman 
et al. [4, 5, 6]. 

It is readily seen that the assertion of Theorem 3.2 follows from Theorem 3.1 with the help 
of Theorem 6.2. 

Summarizing, one can say that the essential equivalence of various forms of Anderson local¬ 
ization (decay of the GFs, EFs, EFCs) is firmly established by now for a large class of random 
Hamiltonians. 

7 Lower regularity of the disorder 

Theorem 7.1. The results of Section 4 remain valid for the marginal probability distributions 
with continuity modulus £y(-) satisfying the following condition: for some C G (0, +oc) and 
an appropriately chosen C > 0, for all e E (0,1/2) 



(7.1) 


Proof Consider first the case where Y k+1 = \_Lj\, hence L k+ 1 > CL]j rT , r > 0. 
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The regularity of the marginal distribution of the random potential V must be sufficient for 
proving a Wegner-type estimate 

p {l|G Bll (£)ll > £?} < 

where /3kSk, replacing 8s k used in the previous section, has to be compatible with our main 
estimates. Denoting e k = L k Sk , we thus should have 

p{||G Bli (E)||> £t }< e ' J *. 

Up to some inessential factors (depending on L k ), the above estimate can be inferred in a stan¬ 
dard way from the continuity of the marginal PDF F v with the continuity modulus of the form 

sy(e) < Ce c 

Next, observe that one has e^ 1 = L s k k < e ClLk : indeed, our estimates by L)~ bk and L k Sk 
are not truly exponential in L k (although that would be very welcome), so we only have e k ~ 

e ±L k lnLk = e ±Lk . Thus 

In In In = In In In L s k k < lnln(ciLfc) < lnln(c 2 I/Q ) < c 3 k. 

At the same time, with (3 k = we have In 8 k 1 > c 4 k, hence one can proceed with the 

scaling algorithm even in the case where 

We conclude that the Wegner-type estimates compatible with the adaptive scaling scheme em¬ 
ployed in Section 4 can be inferred from the following condition upon the continuity modulus 

s y : 

Q 

Sy( £ ) < C £ ln I lne i, 

which is - just marginally - weaker than Holder regularity of any positive order. Pictorially, it 
can be qualified as Holder continuity of ’’almost zero” order. 

The proof in the general case can be reduced to the above analysis, since the double¬ 
exponential growth L k ~ L'l takes over the exponential one, L k = L 4 Y(‘ , after a finite number 
of steps if = if(/i 0 ). Observe that all intermediate calculations and bounds can be re-written 
in terms of strict inequalities (for this is the case with the principal hypothesis, p 0 < a) 2 ), and 
these strict inequalities can be preserved by replacing 8 = Const with 8k = Const/(I + n) k 
during the if steps, provided k > 0 is small enough - depending of course on it. The auxiliary 
constants clearly depend upon the proximity of po to the Germinet-Klein threshold 841 -d . After 
it steps, one can start the scaling procedure with L' 0 := L&. In fact, this would be very close in 
spirit to the Germinet-Klein first bootstrapping step. □ 
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Appendix A Proof of Lemma 1 


Remark 2. The reader may find the proof of the bound (3.5) given below excessively technical 
and opaque. Indeed, we aim here to minimize the size of the ’’double-buffer” nonsingular zone 
around singular cubes; this ’’buffering” technique essentially goes back to [19]. However, a 
simple geometrical consideration shows that, at least for some generous buffering size, resulting 
in some large constant, say, A = 100, one can achieve a non-optimal analog of (3.5) with the 
RHS of the form L k hk< ' Yk+1 ASh ^ ] [ pp; s wou [d still be sufficient for the existence of some scale- 
free probabilistic threshold po G (0, +oo ) for the onset of localization. Below some efforts are 
made to keep p 0 < 29~ M , as the latter appears in the Germinet-Klein BMSA scheme from [26]. 

Fix k > 0 and consider the cube B = B Lk+1 (u), where L k+ 1 = Y k+1 L k , Y k+X = 2K k+1 + 1, 
L k = 34, and its 4-skeleton graph B. For r > 0, let B r = B r {u) be the balls {c G C k : 
d C k{u,c) < r}, and C r := {c G B : ds(ii, c) = r} its boundary (a spherical layer). Then 
B = Br with R = 3K k+ i + 1. Recall that the vertices of B represent the 4-cells (or their 
centers) in the original lattice 7L d . 

We will reduce our analysis of the function B 3 c vp |G B (c, y, E)\ to that of a monotone 
function of one integer variable 

F : r t-3 max |G B (c, y\E)\; r G / := [0, R] . (A.l) 

c:dg(w,c)<r 

We shall see that of particular interest is the sub-interval 

I := [ Kk+h Kk +i + n +1 — 2] = [R — Y k+ 1 , R — 2], 

with R = 3 K k + 1. The reason is that this interval of radii corresponds to the shell of the cube 
B/, :+ i (u) through which the GFs propagate, from the core to the boundary, and it is the decay 
across this spherical layer that we have to assess. This analysis is essentially reduced to the 
points with r > R — Y k+ p. for smaller radii, we apply first the GRI and jump to the boundary of 
the core (with maximal penalty upper-bounded thanks to the CNR condition), and then start the 
analysis presented below (cf. item (C)). 

Introduce another function, 


/ : r max \G B (c,y, E)\ = max |G B (c, y E)\ , 

c:dg(u,c)=r 


(A.2) 


so that F(r) = max r /< r f(r'). 

Call a vertex c G B non-singular if the associated ball B Lk (c) C Z d is (E,Lf bk )- NS, and 
singular, otherwise. Respectively, call r G / non-singular if all vertices c with dist^w, c) = r 
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are non-singular, and singular, otherwise. The notions of singularity/non-singularity do not 
apply to r 6 [R — 1, R — 2]. 

By definition of the non-singular balls (cf. Definition 2 and (2.11)) we have the following 
inequalities: 

(A) for any non-singular r G I, 


f(r)< Cw,k L k 

re[r-l,r+l] 


/(O; 


(B) owing to the assumed CNR-property of B ifc+1 (u), for any r < r' < R — 2, an application of 
the GRI (cf. (2.17)) to the ball A r 'L k (u) C B Lk+1 (u ) gives 


/(r)<G^ + 1 /(r'). 


(A.3) 


To be more precise, an application of the GRI is required for r < r' — 1, while for r — r' the 
inequality (A. 3) follows trivially from Y^ +1 L 8 k k +1 > 1. 

Combining (A) and (B), we come to the following statement: 

(C) Assume that for some r' G /, all points p G [r' + 3, r' + 5] are non-singular. Then for all 

r G [0, r' + 5] 

Hr) < Cp k F(r' + 6 ), (A.4) 

Notice that for r = r' + 5, (A.4) follows immediately from the assumed non-singularity of 
the point r' + 5, so it remains to be established only for r < r' + 4. 

For the proof, we first apply (B): 


F{r' + 4) = max f(p) < C w , k L s * x f(r' + 4). (A.5) 

p<r'+ 4 

Next, apply (A) to r' + 4 (which is non-singular by assumption): 


f(r' + 4)<C^ k Lp- 


max f(r"), 
r"e[r'+3,r'+5] 


(A.6) 


thus 


F{r' + A)<C w \L s * +l L k 


-bk 


max f(r 
' €\r'+3^'+5\ 


n\ 


(A.7) 


Apply (A) once again to the three points r" G [r' + 3, r' + 5] (all of which are non-singular by 
assumption): 


max f(r")<C Wk .L . k max max f{F") 

r"e[r'+3,r'+5] ’ r"G[r'+3,r'+5] r'"6[r"—l,r"+l] 


^ /o—1 r —b k 

A 

if—1 T ~bk 


max /(r w ) 


r /// G[r / +2,r / +6] 


F(r' + 6 ). 


(A.8) 
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Figure 5: In this example, K = 7, Y = 2 • K + 1 = 15, S = 1, N = Y — 5S — 1 = 9, R = 22. It is assumed 
that there is no pair of disjoint singular intervals of the form [p — 1 , p + 1 ], r% is the smallest integer in / which 
is singular; it is the radial projection of the center c of a singular ball in the skeleton graph. It is this minimality 
property which implies that r 2 = — 1 must be non-singular, despite the fact that the intervals [r^ — 1 , r 2 + 1 ] and 

[r 3 — 1, r 3 + 1] overlap. On the other hand, due to the overlap of [r 3 — 1, r 3 + 1] with [r 3 , + 2], the point r 3 + 1 

may (or might) be singular, without producing a disjoint singular pair. Therefore, we still can use the property (A) 
starting off the point (and aiming at r 2 + 1 = r 3 ), but leaving from ?~ 3 , we have to make a longer flight with 
possible ’’destinations” (i.e., reference points) ranging in [(V 3 + 4) — 2, ( 7-3 + 4) + 2] =[r 3 + 2, 7-3 + 6], The longest 
flight consumes the distance 6, instead of 1 that we would have for a non-singular departure point; this results in a 
loss of 5 points. The thick gray intervals indicate the ’’good” points providing the factors q < (3F — l)~ d L~ b < 1 
in the ’’radial descent” recursion: F(ri- 1 ) < qF(ji). The point rg is used as the last reference point, but we can 
only bound F(rg) by the global maximum of F, since the GRI cannot be applied at a center c of the skeleton graph 
B with d b(u,c) > R — 1. In this example, we have the guaranteed decay bound F(rg) < q 9 F(rg) < q 9 F(R), 
with 9 = Y - 55 - 1. 


Collecting (A.7) and (A.8), the assertion (C) follows, since F(r) < F(r') for r < r'. 

Now pick any maximal collection of disjoint singular cubes B L k (ci), i = 1,... ,n < S, 
denote pi = d B (u,Cj), and associate with each c, an interval [5,, /3J = [p, — 1, p t + 5]. Next, 
decompose the union of intervals [5j, 13i\ into a disjoint union of maximal non-overlapping in¬ 
tervals Ji = [oti,/3i\, 1 < i < n' < 7i, so that fa < a i+ 1 — 1; the equality a i+ 1 = fa + 1 is 
permitted. 

Note that for any i, all points r e [/3i — 3, /3i\ are non-singular, otherwise we would have to 
augment J, t by fusing it with the interval [r — 1, r + 5] overlapping with J t , which contradicts 
the maximality of Ji. 

Let V = I\UiJi and enumerate the points of/'in the natural increasing order: r 0 ,r , 1 ,.... Tm- 
Next, add formally the last point tm+ 1 R — 1: it will be used as the last reference (desti¬ 

nation) point, although it cannot serve itself as a departure point. If r t is non-singular, then we 
have r i+ \ = r* + 1 and 

Jdn) < C^ k Ll b -F(r w y, 

otherwise, we can apply (C) and obtain 

F(r<) < = \\W\\~%7 d L?"* Ll^F^) 

< HlVir 1 (L~ k b -+‘-Y‘l ,) ■ Y k - + \ L~ t b -F( r i+I ), <A ' 9> 
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Figure 6 : Here K = 7, Y = 2 • K + 1 = 15, S = 2, N = Y - 5S - 1 = 4. Three variants of the 
worst case scenario, with non-overlapping intervals J \, J 2 . The ’’useful” points of I, serving as departure points 
and providing a small factor are drawn as thick intervals (black, if r* is singular, and dark gray, otherwise) and 
numbered as ro,..., rs- Here we have the guaranteed decay bound F(ro) < q A F{r,{) < q A F(R). One can see that 
the number of useful points becomes larger if two or more intervals Ji overlap with each other, or if at least one of 
them overlaps with the complement of / = [R — Y, R — 2] ( = [7, 20], in this example). 


For k = 0, we have b 0 — s 0 = s 0 — 4 = f3 G (0,1], while Y 1 — 9 < L 0 , so 


P P' 

n \ -It f±H 

L ~ {b °~ S0) Yl S0 = Lq P Yi ^ < 1 , 

1 +;' 


owing to the assumption L 0 > Yl v (cf. (3.2)). For k > 1, we have s k = | b k , and again, 

T —b k +S kv Sk _ T~ 6 Sk V^k tz' T~ 6 Sfc+ 8 Sfc ^ 1 
X k+1 — 1 k+1 — L i 

hence for any k > 0 , 


F( ri ) < WWW^L^Fir^) . (A.10) 

Collapsing the intervals J % removes from / at most 5S k+ 1 points (see Fig. 6 ). This upper bound 
becomes sharp if the radial projections of all singular /./,-balls in the collection (fixed at the 
beginning) are non-overlapping. Hence \I'\ > (Y k+] — 1) — 5,S'/,. + , = N k+1 , so we obtain 


m < c w \l-j n *+' < 


— 1 ‘-‘k 


1 V d 

L f—bkJk+l ^ 1 fc+1 T~ b kJk+l 

V d L k - r L k 

1 k Fw,k +1 

— L k ^w,k+ 1 


(A.ll) 


—bkNk+i ^ j-d/8 l r — frfc-N/e+i 

J-Jl, OlI7L I 1 L 


We conclude that 


Lemma 1 is proved. 


Gb, f(0) < if if "** 1 ■ 


(A. 12) 
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